Dimensional Reduction in Bose-Condensed Alkali-Metal Vapors 
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We investigate the eflects of dimensional reduction in atomic Bose-Einstein condensates (BECs) 
induced by a strong harmonic confinement in the cylindric radial direction or in the cylindric axial 
direction. The former case corresponds to a transition from 3D to ID in cigar-shaped BECs, while 
the latter case corresponds to a transition from 3D to 2D in disc-shaped BECs. We analyze the first 
sound velocity in axially-homogeneous cigar-shaped BECs and in radially-homogeneous disc-shaped 
BECs. We consider also the dimensional reduction in a BEC confined by a harmonic potential both 
in the radial direction and in the axial direction. By using a variational approach, we calculate 
monopole and quadrupole collective oscillations of the BEC. We find that the frequencies of these 
collective oscillations are related to the dimensionality and to the repulsive or attractive inter-atomic 
interaction. 



PACS Numbers: 03.75.Kk 

In the experiments with Bose-Einstein condensates 
(BECs) the atomic cloud is confined by a magnetic or 
magneto-optical trap. In the simplest case the confin- 
ing field acts as an anisotropic harmonic potential. Re- 
cent experiments [1,2] have shown that it is possible to 
change the frequencies of the harmonic confinement in- 
ducing a dimensional reduction of the system: from 3D 
to ID and from 3D to 2D. Effects of reduced dimension 
in BECs have been theoretically investigated by many 
authors [3-6], but less attention has been paid to the 
transition from a dimensional regime to another. 

In this brief report we study in detail the di- 
mensional reduction in BECs by using nonpolynomial 
Schrodinger equations, which are derived from the 3D 
Gross-Pitaevskii equation (CPE) [7]. In this way we 
calculate chemical potential, elementary excitations and 
first sound velocity for cigar-shaped and disc-shaped 
BECs. Moreover, we apply a variational technique [8] to 
obtain monopole and quadrupole collective excitations 
of the BEC with both positive and negative scatter- 
ing length as- In particular, we analyze the effects of 
anisotropy and inter-atomic strength on these collective 
frequencies. 

From 3D to ID. - We have recently shown that the dy- 
namics of a cigar-shaped BEC of N atoms, confined by 
a harmonic potential of frequency uor in the cylindric ra- 
dial direction (x, y) and by a generic potential V{z) in the 
cylindric axial direction z, is well described by an effec- 
tive ID nonpolynomial Schrodinger equation (ID NPSE) 
for the axial wave function /(z, t) of the BEC [7]. Within 
the ID NPSE approach, the BEC is approximated by a 
Gaussian in the cylindric radial direction. The Gaus- 
sian width a is related to / by ct = (1 + gA^I/P)^/'', 
where g — 2as/ar is the effective interaction strength 
with Or = h / (muJr) the characteristic harmonic length 
in the radial direction [7]. Note that we write length is 
in units and energy in units hujr- 



In the limit of weak radial confinement, which corre- 
sponds to the 3D regime, the axial density satisfies 
the condition gN\f\'^ >> 1 and the ID NPSE reduces to 
a ID nonlinear Schrodinger equation with quadratic non- 
linearity. In the opposite limit gN\f\'^ << 1 of strong 
radial confinement, which corresponds to a ID regime, 
the ID NPSE becomes a ID nonlinear Schrodinger equa- 
tion with cubic nonlinearity, also called ID CPE [4]. It 
is important to remember that the ID NPSE is derived 
from the 3D GPE, but the 3D CPE (and also ID NPSE) 
is no more valid in the limit of strong diluteness: in this 
case the atomic cloud becomes a Tonks gas, i.e. a gas of 
impenetrable Bosons [9]. In our units the system is far 
from the Tonks regime if the condition (?iV|/p ^ g^ is 
satisfied [10]. 

In previous papers [7] we have verified that the ID 
NPSE is able to capture the full 3D regime, the full ID 
regime (apart the Tonks gas) and the region in between. 
Here we show that an analytical expression for the shape 
of the axial density can be obtained by neglecting the ki- 
netic term, the so-called Thomas-Fermi (TF) approxima- 
tion, in the ID NPSE. Then, setting f{z,t) = e"^''*/o(z) 
one finds for the axial density profile pi{z) = N\fQ{z)\'^ 
the following expression 



pi{z) = — yn{zf - 3 + fi{z)y^fi{zy^ + 3 



(1) 



where /i(z) — fi — V{z) with fi the chemical potential. 
This density profile reduces to pi{z) — (2/i(z))^/(9g) 
when gpi >> 1, that is very close (with the factor 4/9 
instead of 1/2) to the axial density profile of the TF 3D 
regime obtained from the 3D GPE. Eq. (1) reduces to 
pi{z) = {p,{z) — 1) / g when gpi << 1, that is exactly the 
density profile of the TF ID regime obtained from ID 
GPE. 

In the case V{z) — 0, which corresponds to a homo- 
geneous BEC in the z direction with periodic boundary 
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conditions, one gets the axial velocity of first sound as 
ci = {pi4^Y^^ form which one finds 



9Pi _ 1 9 Pi 



'4^1 + gpi 



4(l+5pi)3/2 



(2) 



In Figure 1 we plot ci as a function of gpi. The figure 
shows that ci strongly depends on the dimensionality. 
In fact, our expression for the velocity of sound gives 
ci = ■\/3/4(5/9i)^/'* in the 3D regime and it becomes 
ci — \fgp\ in the ID regime. Note that our ID value 
for ci is in full agreement with the result of Menotti and 
Stringari [10], while at higher density our result c\ differs 
from the exact one obtained by the 3D GPE [10], only for 
the constant that is v^/S instead of 3/4. This small dis- 
crepancy (relative error less than 3%) in the 3D regime is 
a consequence of using in the ID NPSE a Gaussian radial 
wave function. Thus, our Eq. (2) for the sound velocity 
Ci interpolates very well between the exact rcsiilts of low 
density (ID regime) and high density (3D regime). 




FIG. 1: First sound axial velocity ci versus intcr-atomic 
strength gp\, where g = 2as/ar with Os scattering length and 
Or radial harmonic length, pi is the axial density of the cigar- 
shaped axially-homogeneous BEG. Length in units ar, density 
in units l/ur and time in units l/wr- 

In the case V{z) = muj1z^/2 the low-frequency ele- 
mentary excitations, which are collective oscillation of 
the BEG, can be more easily calculated with a varia- 
tional approach [8] from the full 3D GPE instead of us- 
ing the ID NPSE. The most natural choice is a Gaussian 
wave function with two variational parameters: the ra- 
dial width cr and the axial width ry. Then the effective 
potential energy U of the atomic cloud obtained from the 
3D Gross-Pitaevskii energy functional is given by 
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where A = u)r/^z is the anisotropy parameter of the 
confining harmonic potential. With A > 1 the BEG is 
cigar-shaped, while for A < 1 the BEG is disc-shaped. 
The equilibrium point {(j*,rf), corresponds to the min- 
imum of the effective potential energy of the system. 



The low-frequency collective excitations of the BEG are 
the small oscillations around the equilibrium point. For 
A = 1 these two normal mode frequencies can be identi- 
fied as the monopole frequency lom and the quadrupole 
frequency ujq [8]. One finds that for A >> 1 the fre- 
quency Wm is mainly transverse and the frequency loq 
is mainly axial, while for A << 1 the frequency lom is 
mainly axial and the frequency uq is mainly transverse. 
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FIG. 2: Monopole frequency ujm and quadrupole frequency 
(jOq of collective oscillations for a cigar-shaped BEC. Case of 
positive scattering length Os. \ = oJt/ijJz is the anisotropy of 
the harmonic trap, uh = ^/fi/{mijJH) with ujh = (w^w^)^^^. 

In Figure 2 we plot um and uiq as a function of 

the inter-atomic strength for different values of the 
anisotropy parameter A. In the spherical case (A = 1) 
the monopole frequency um ranges between the value 
(lomI^t)^ = 5 of the 3D regime of mean-field bosons 
(isotropic TF 3D) and the value [ujM/'jJrY = 4 of quasi- 
ideal bosons. By increasing A the frequency lom becomes 
a constant, given by {ujm /ij^r)"^ = 4, independently on 
the interaction strength. 

The transition from 3D to ID is even better shown 
by the quadrupole frequency loq (sec also [12]). In this 
case, for very large values of A the frequency ujq shows 
three plateaus: at [ujq/uJz)'^ = 5/2 the BEG is in the 3D 
regime (anisotropic TF 3D), at {ujq/uJz)'^ = 3 the BEG is 
in the regime of ID confinement and mean-field bosons 
(TF ID), and at {ljq/u;^)^ = 4 the BEG is in the regime 
of ID confinement and quasi-ideal bosons (free bosons 
ID). Thus Figure 2 shows that it is necessary a very 
strong anisotropy A to reach the regime of ID confine- 
ment and mean-field bosons, and ujq is able to discrim- 
inate between this regime and that of ID confinement 
and quasi- ideal bosons. Note that, as shown in [10], the 
Tonks regime is avoided if N(as/aH) ^ A^/'^(as/a/f)'^: 
the results of Fig. 2 are reliable for all values of horizon- 
tal axis under the condition ttg/ajj <^ A~^/^10~^'^/^. 

Our variational approach can be applied also to the 
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case of an attractive BEC (ag < 0) . In this case it exists 
a critical inter-atomic strength below which there is the 
collapse of the BEC [11], as shown in Figure 3. 
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FIG. 3: Monopole frequency ujm and quadrupole frequency 
ujQ of collective oscillations for a cigar-shaped BEC. Case of 
negative scattering length as- 

In the spherical case the monopole frequency lom grows 
by decreasing the interaction strength up to the collapse, 
but it is sufficient a small anisotropy to freeze com to a 
constant value {ujM/'^rf' = 4. Instead the quadrupole 
frequency loq strongly increases with A but always goes 
to zero at the collapse. 

From 3D to 2D. - The dynamics of a disc-shaped BEC 
of N atoms, confined by a harmonic potential of fre- 
quency Wz in the cylindric axial direction z and by a 
generic potential W{x,y) in the cylindric radial direc- 
tion is well described by an effective 2D non- 
polynomial Schrodinger equation (2D NPSE) for the ra- 
dial wave function (j){x, y, t) of the BEC. Within the 2D 
NPSE approach, the BEC is approximated by a Gaus- 
sian of width rj in the cylindric axial direction, given by 
rf — ')N\(l)^r} —1 = 0, where 7 = 2\/27ras/az is the effec- 
tive interaction strength and = \Jfij (muiz) the char- 
acteristic harmonic length in the radial direction. Here 
length is in units and energy in units fuVz- 

Note that in the limit of weak axial confinement, 
which corresponds to the 3D regime, the radial den- 
sity N\(l)\^ satisfies the condition 7A''|<j!)|^ >> 1 and the 
2D NPSE reduces to a 2D nonlinear Schrodinger equa- 
tion with quadratic nonlinearity. In the opposite limit 
jN\(l)\'^ << 1 of strong axial confinement, which corre- 
sponds to the 2D regime, the 2D NPSE becomes a 2D 
nonlinear Schrodinger equation with cubic nonlinearity, 
also called 2D GPE [7]. 

In the case W{x,y) = 0, which corresponds to a homo- 
geneous BEC in the cylindrical radial direction with pe- 
riodic boundary conditions, the chemical potential reads 



It is easy to show that Eq. (4) reduces to /x = 
(5/4) in the 3D regime (77 = (7^2)^''"^) and it be- 

comes /i = gp2 + 1/2 in the ID regime (ry = 1). The 
term 1/2 is the axial energy in scaled units. From the 
2D NPSE one calculates the radial velocity C2 of first 
sound as C2 = {pi^^f^- 
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FIG. 4: First sound radial velocity C2 versus inter-atomic 
strength 7^2, where 7 = 2\/2itas/az with as scattering length 
and az radial harmonic length. p2 is the radial density of the 
disc-shaped radially-homogeneous BEC. Length in units az, 
density in units 1/a^ and time in units 1/oJz- 

In Figure 4 we plot C2 as a function of 7^2- The velocity 
of sound is C2 = (5/6)(7p2)^/^ in the 3D regime and it 
becomes C2 = ^7^2 in the 2D regime. Our 2D value 
for C2 is in full agreement with the result obtained from 
2D GPE in the TF approximation, while the 3D value 
of C2 differs from the exact one, obtained from the 3D 
GPE, only for the constant that is Z'^I^-k^I'^ jiJ 1^ instead 
of 5/6. This error (less than 2%) in the 3D regime is a 
consequence of using in the 2D NPSE a Gaussian axial 
wave function. Thus, our sound velocity C2 interpolates 
remarkably well between the exact results of low density 
(2D regime) and high density (3D regime). 

In the case W{x,y^ = muj^{x^ + y^)/2 the collective 
oscillations of the BEC can be obtained from Eq. (3) 
with the condition of disc-shaped BEC, i.e. < A < 1. 
In Figure 5 we plot monopole and quadrupole frequen- 
cies of the; BEC with a weak anisotropy. The monopole 
frequency lom does not depend on A in the limit of very 
small inter-atomic strength and coincides with the ideal 
gas limiting value. In the regime of mean-field bosons, 
wm changes with A from {wm/i^z)'^ = 5 (isotropic TF 3D) 
to (ujm/^z)^ = 3 (anisotropic TF 3D). In a similar way 
the quadrupole frequency loq moves from (wq/oJj.)^ ~ 2 
(isotropic TF 3D) to {ujQ/cOrf = 10/3 (anisotropic TF 
3D), as predicted by Stringari [12]. 
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FIG. 5: Monopole frequency com and quadrupolc frequency 
ujQ of collective oscillations for a disc-shaped BEG. Gase of 
positive scattering length Os with small anisotropy. 

By further increasing the disc-shaped anisotropy the 
main effect is that of shifting the transition from the ideal 
gas limit to the 3D regime at higher values of the inter- 
action strength. As shown in Figure 6, both monopole 
and quadrupole frequencies display this behavior. Thus, 
these frequencies do not discriminate the regime of 2D 
confinement and mean-field bosons from the regime of 
2D confinement and quasi-ideal bosons. 
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FIG. 6: Monopole frequency ujm and quadrupole frequency 
uiQ of collective oscillations for a disc-shaped BEG. Gase of 
positive scattering length Os with large anisotropy. 

The frequencies simply give a signature of the transi- 
tion from ideal gas to TF 3D. Clearly, under strong con- 
finement it is necessary to increase the inter-atomic inter- 
action to mantain the 3D regime. Note that for the TF 
2D regime the dispersion relation of collective frequen- 
cies can be analytically obtained using the same hydro- 
dynamical approach applied by Stringari for the TF ID 
regime [12]. In this way one recovers the result = Au)"^ 
for the quadrupole frequency. 

The calculations with negative scattering length ag 



suggest that monopole and quadrupole oscillations of a 
highly deformed disc-shaped BEC do not depend on the 
inter-atomic strength. 
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FIG. 7: Monopole frequency ujm and quadrupole frequency 
uiQ of collective oscillations for a disc-shaped BEG. Gase of 
negative scattering length Os. 

As shown in Figure 7, only very close to the collapse 

ujQ goes to zero while lum remains essentially constant. 

In conclusion, our results exactly reproduces the ones 
obtained in specific dimensional regimes with the hydro- 
dynamic method or the sum-rule method, but, they also 
give the crossover from one regime to another. 
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